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Abstract 



Statistical properties of optical fields in nonlinear couplers composed of two waveguides 
in which Raman or Brillouin processes (with classical pumping) are in operation and which 
are mutually connected through the mutual Stokes and/or anti-Stokes linear interactions are 
investigated in the framework of a generalized superposition of coherent fields and quantum 
noise. Heisenberg equations describing the couplers are solved both analytically under special 
conditions and numerically in general cases. Regimes for nonclassical properties of optical 
. fields, such as sub-Poissonian photon number statistics, negative reduced moments of inte- 

grated intensity, and squeezing of quadrature fluctuations are discussed in cases of single and 
compound fields. General results are compared with those from the short-length approxima- 
Q\ • tion. 

1 Introduction 

Recently increasing attention has been devoted to the study of nonlinear couplers composed of 
two or more waveguides with mutually connected modes by means of evanescent waves [jlj, @, [5| . 
Also quantum statistical properties of such devices have been studied j|, ^ . Especially, couplers 
composed of linear and nonlinear waveguides have been examined for both directional ^| and 
contradirectional |(| Q propagation of light beams. Also couplers composed of two nonlinear 
waveguides have been investigated || . Various nonlinear optical processes were assumed to operate 
in nonlinear waveguides, such as second harmonic generation S (and references therein), frequency 



downconversion (|, |To[ [n| [l2j and Kerr effect p], [l4| p| 

We examine quantum statistical properties of couplers composed of two waveguides in which 
Raman or Brillouin processes are active and in which modes (Stokes and anti-Stokes) mutually 
linearly interact through evanescent waves. We suppose classical strong laser pumping of the 
nonlinear processes in both the waveguides. Then we can apply the method of the generalized 
superposition of coherent fields and quantum noise to study spatial dependence of statistical quan- 
tities, such as photon number distribution, reduced moments of integrated intensity, quadrature 
variances and principal squeeze variance for both single and compound modes. Special attention is 
devoted to nonclassical properties of light exhibited by sub-Poissonian photon statistics, negative 
reduced factorial moments or squeezing of quadrature variances below the coherent-state values. 
Statistical properties of modes involved in Raman or Brillouin processes have been studied for a 
long time (for reviews, see, e.g. 16, |l7|, Q). Especially, results contained in papers |o[ ^lf , 



in which quantum statistical properties of Brillouin ||19|, |20| and Raman pll processes were inves- 
tigated in the framework of the generalized superposition of coherent fields and quantum noise, 
can be employed. The investigated model represents a generalization of the above models in the 
direction of incorporation of the influence of linear Stokes and anti-Stokes interactions. 
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Figure 1: Scheme of interactions in the quantum nonlinear coupler formed from two nonlinear 
waveguides operating by means of Brillouin or Raman processes, in which Stokes and anti-Stokes 
modes mutually interact through evanescent waves; Sj, Aj, Lj, and Vj denote Stokes, anti-Stokes, 
laser pump, and phonon modes, gs } (g A j ) is Stokes (anti-Stokes) nonlinear coupling constant in 
the jth waveguide, and k$ {^a) represents Stokes (anti-Stokes) linear coupling constant. 



Application of some earlier results for Raman scattering to obtain sum and difference squeezing 



of vacuum fluctuations were given in ]22| , 23 , including the quantum correlations between the 
Stokes and anti-Stokes fields, whereas the correlations between the pump and Stokes waves were 
considered and measured in ]24| ]. Raman effects can also play important role in quantum soliton 
propagation p5[ | and in solving various inverse scattering problems [ p6[ . 

The dynamics of the couplers is solved analytically, numerically, and in a short-length ap- 
proximation in Section 2. Quantum statistical properties of modes of the couplers are described 
in Section 3. Section 4 contains a discussion of the statistical behaviour of the systems under 
discussion based on results of the numerical analysis. Conclusions are summarized in Section 5. 



2 Dynamics of nonlinear couplers 



A model of the coupler composed of two mutually interacting waveguides with Raman or Brillouin 
processes is quantally described by the following momentum operator (for a scheme of involved 
interactions, see Fig. 1): 



G = 



J~] hkjalcbj + hg Al a Ll a Vl a Ai +hg Sl a Ll alr i al i +h.c. 

j=Li,Si,AuVi 



j=L 2 ,S 2 ,A 2 ,V2 



t At 



hg A2 a L2 a V2 a\ + hg S2 a L2 a' v a 



V 2 U, S 2 



[hUsas^^ + TihAa,A l a A2 + h.c. 



(1) 



where aj (a]) are annihilation (creation) operators of laser (L), Stokes (S), anti-Stokes (A), and 
vibration (V) modes in the first (1) and the second (2) waveguides. Wave vectors of the corre- 
sponding modes are denoted as kj. Coefficients gg 1 and gs 2 {gAi an d 9a 2 ) describe nonlinear 
Stokes (anti-Stokes) interactions in the waveguides 1 and 2. Coefficients k$ an d k A are related to 
linear coupling between Stokes modes and anti-Stokes modes, respectively, in both the waveguides 
by means of the evanescent waves. The symbol fi denotes the reduced Planck constant and h.c. 
means Hermitian conjugate terms. 

The momentum operator G in (1) reflets a symmetry between the two waveguides, which is 
based on the exchange of the suffices (Li, Si, A\, V{) and the suffices (L 2 , S^, A 2 , V2) with the 
simultaneous exchange of constants ks, ka with their complex conjugate ones k* s , k* A . This 
symmetry has to be conserved in spatial development of operators and all quantum statistical 
quantities, so we restrict our discussion of results only to modes in waveguide 1; results valid for 
waveguide 2 are obtained simply from the above mentioned symmetry. 
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Quantum spatial development of an operator a is described by the Heisenberg equation ih^ = 
[G,a], where z is a spatial variable along the propagation direction, G is the above introduced 
momentum operator and the symbol [ , ] denotes a commutator. The Heisenberg equations for 
the operators Aj in the interaction picture [Aj = aj cxp(— ikjz)] and under the assumption of 
strong classical laser pump modes L± and L 2 (aj — ► aj exp(ifcjz) for j = L\,L,2, aj being classical 
amplitudes) have the form 



dA Sl 
dz 

dA Al 
dz 

dA Vl 
dz 



igStA^ +in* s A S2 , 
ig Al A Vl + in* A A A2 , 



(2) 



where 



9S 3 

K A 



9s j Q.L j exp(iAk S] z), 
g A] a Lj exp(iAk Aj z), 
k s exp(iAK s z), 
k A exp(iAK A z) 



for j = 1,2, 



(3) 



with the phase mismatch vectors Akg^ = ki Jj — kv d — ks j , Ak Aj — fe^, + ky 6 — k Aj for j = 1,2 and 
AK S = k Sl - k S2 , AK A = k M - k A2 . 

The Heisenberg equations (2) are linear, which enables us to solve them simply numerically 
and to retain operator character of the solution (see the next section). They can be also solved 
analytically under special conditions. Furthermore a short-length solution can be reached. These 
three kinds of solutions will be described in the following subsections. 

The relation 



d_ 

dz 



(4, (z) + A Aj (z)A Aj (z) - Al (z)A Sj (z)) 



i=i,2 



= 0, 



(4) 



which represents a conservation law for photon numbers, can be derived directly from (2). 
2.1 Numerical solution 

The set of equations (2) can be conveniently written in the matrix form reflecting the symmetry 
of the model: 

rfA -™ a 
— = iMA, 

dz 

(5) 

where the vectors Ai, A 2 and the matrices Mi, M 2 , Mi 2 are defined as follows 

°J 

it 



A = 


' Ax " 


,M = 


Mi M12 " 




. A 2 _ 




. M* 12 M 2 _ 



Ai = 



A Aj 

A Vj 
A f 
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for j = 1,2, 



(6) 



The above set of equations is composed of two independent sets of equations for operators 
(A' Si , Aa x , Av ± , Ag , Aa 2 , Ay 2 ) and for their Hermitian conjugates. It is more convenient to ma- 
nipulate with only one larger set of equations, because the procedure to solve them can be easily 
generalized for the case of losses in modes, in which also fluctuating Langevin forces must be 
taken into account; this results in coupling of operators and their conjugates through correlation 
functions of fluctuating forces (see, e.g. p7|). 

Solution of equations (5) can be written in the form A(z) = exp(iM2:)A(0) supposing that the 
matrix M is z— independent. Since we assume in the following that all mismatches are zero, the 
matrix M is really z— independent. Diagonalization of the matrix M then enables us to determine 
spatial dependence of the evolution matrices U and V in the operator solution (for details, see e.g. 
@) 

6 

Aj(z) = ]T [u jk (z)A k (0) + V jk (z)A{(0)\ for j = 1, . . . , 6, 

fc=i 

where the indices j, k run over the symbols S±, A\, V\, S2, A2, and V2. 



(7) 



2.2 Analytical solution 

In a special case, characterized by the conditions 

IffSil = \9s 2 \, \9Ax\ = \9a 2 \, M = \ka\, 
k\ 9a 2 _ n s g* S2 



\ka\9A! \Ks\g*s 1 , 

we can solve analytically the set of equations (5) using the following transformations (k = \ks\) 



(8) 



bs - 

cs = 

b A - 

CA = 



cv 



—j= exp(— inz) 

— = exp(iKz) 
V2 

— = exp(— inz) 

—= expfiKz) 
s/2 



1 

1 

71 



exp(— ikz) 
exp(iKz) 



A Sl 
A Sx - 

A Al 

A Al - 

A Vl 



__S_ A 

I I ■ 



Si 



^A 



s 2 



A A f 



A 



Vi 



\K A \ 
^ K A 9A 2 

I«a| 9a 1 

K* A 9A 



A 



v 2 



\ka\ 9a 1 



A 



v 3 



(9) 
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The application of the transformations (9) results in the following sets of equations: 



d_ 

dz 

d_ 

dz 



' b A ' 










9 At 




r b A ' 


C S 


= i 












c s 


. by . 






9Si 


— K 




. by _ 


CA 














CA 




= i 








-9* Sl 
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.9M 


ffSl 


K 







(10) 



and two sets obtained by Hcrmitian conjugation of the above equations. These equations can be 
solved analytically (e.g. by finding eigenvalues and eigenvectors) and the inverse transformation 
leads to the solution of the operator equations (5) in the form: 



4W = 



The abbreviations 



1 



+ 



|ffAi| 2 chh 

ff^ffSi 



Iff Si 



)]4(0) 



chh - — sh si* + ch cl* 

I* 



+ i 



+ i 



1J± 

r 2 \ks\ 



■9a 1 9s 1 k a 



\ K A\ 



j^-^shsf + chcr 
"4,(0) 

5Al | 2 shh+|. 9Sl | 2 (^ch sf+sh cl*) 

4 2 (0) 



4(0) 



shh — — ch si* — sh cl* 

I* 



+ 2z^ch sl*it (0) - 2^^sh sl*it (0), 



I* ~ VlV ' | Kg | I* " ' " V ~ 



AaA*) = 



9A 1 gS 1 



chh + — sh si — ch cl 



^s 1 (0) + 2*^ch si A Vl (0) 



+ i 



| ff5l | 2 chh+| ffAl | 2 (ysh sl-chcl)] A Al (0) 

\g Sl | 2 shh - \g Al | 2 (ych si + sh cl)] A A , (0) 



— i- 



1^1 

r 2 \k a \ 
■9A 1 gs 1 K S 



r \ks 

K A 9A 2 



shh — -ch si — sh cl 



4,(o) 



\k a \ I 



sh si Ay 2 (0), 



A Vl (z) = 2i^-chslA Si (0) + 2i^-chsll Al (0) 



9a, 



-sh si + ch cl 



4(o) 



+ 2 ^l^shsl4 (0)-2^^ 7 shsli A2 (0) 



I \k s \ 

K* A 9A 2 
\ka\ 9 A x 



— ch si — sh cl 



shh = sin(/cz), 
sh = sin(Kz/2), 
si = sin(lz/2), 

r = VlffsJ 2 - 



.9aJ 2 , 



I \ka\ 

4(o). 



chh = cos(k2), 
ch = cos(k2;/2), 
cl = cos(Zz/2), 
/ = \J k 2 - 4r 2 



(11) 



(12) 



have been used in the above expressions. Expressions for operators A S2 , A^ 2 , and Ay 2 are obtained 
from the expresions (11) using the before mentioned symmetry of the coupler. Spatial develop- 
ment of creation operators follows from Hcrmitian conjugation of the expressions for annihilation 
operators. 
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2.3 Short-length solution 

A short-length solution of equations (5) valid up to z 2 has the form: 
A Sl (z) « A gi +i(g 8l A^+K* s As a )z 

+ (?Si9A x A ] Ai + \g Sl \ 2 A Sl - K* s g S2 Al 2 - \n s \ 2 A Sl ) z 2 /2 
A Al (z) w A Al + i (g Al A Vl + K* A A A2 ) z 

- (sSiSAiA^ + IffAj 2 ^! + K* A g A2 Av 2 + \k a \ 2 A Ai ^ z 2 /2 
A Vl (z) « i Fl + i (g* Al A Al + g Sl A f Si ) z 

- ({\9aA 2 - \gsA 2 )A Vl + g A y A A A2 - .g Sl ^4 2 ) * 2 A (13) 

Operators on the right hand sides of expressions in (13) are taken at z = 0. Expressions for the 
other operators are reached from the symmetry and by means of Hermitian conjugation of the 
above expressions. 

3 Quantum statistical properties of nonlinear couplers 

Quantum statistical properties of optical fields involved in the above processes can be studied in 
the framework of the generalized superposition of coherent fields and quantum noise ( |p"8f , Sees. 
8.5, 9.3, and 9.4), which is characterized by the normal characteristic function CV({/3j}, z) in the 
Gaussian form: 

C U ({0j},z) = ex p|E ~B j (z)\0 j \ 2 + (~Ci(*)# 2 +cc.) 

6 6 

+ E E [D jk (z)(3*p* k + D jk (z)f3j(3* k + o.c.} 

j=l k=l,j<k 

+ pm*(z) -c.c.}^, (14) 

where c.c. denotes the complex conjugate terms. The complex amplitudes £i(z), . . ., £e(z) corre- 
spond to the annihilation operators As 1 {z), A Al (z), Av 1 (z), As 2 {z), A A2 (z), and Av 2 (z), respec- 
tively. The quantum noise functions Bj(z), Cj[z), Dj k {z), and Dj k {z) defined by the relations 

Bj(z) = (AA](z)AA 3 (z)), 

C 3 {z) = {(AA 3 {z)) 2 ), 

D jk (z) = (AA j (z)AA k (z)), 

D jk (z) = -{AA]{z)AA k {z)), (15) 

can be determined in terms of values of these quantities for incident fields, Bj, Cj (Bj = Bj(0) + 1, 
Cj = Cj(0) if antinormal ordering of field operators is adopted), provided that the operator solution 
of the corresponding Heisenberg equations is known (for details, see e.g. p§|]). Antinormal ordering 
adopted for incident field operators enables us to describe nonclassical states of light. Values of 
quantities Bj and Cj are given in the case of the squeezed fields with additional noise by 

Bj = cosh (fj ) + (n chj ), 

Cj = ie X p(^) S inh(2r j ), (16) 



(i 



where rj are the squeeze parameters of the incident beams, 9j are the squeeze phases and (n c hj) 
represent the mean number of external noise particles in the jth mode. In general signal squeezed 
light with an additional noise is described by the above expressions; when rj = 0, (n c hj) = 
(Bj = 1, C } ; = 0) we have a coherent state and the superposition of the signal coherent field and 
noise is characterized by r.j — 0, (n c hj) 7^ (Bj = (n c hj) + 1, Cj = 0). 

The photon number distribution p(nj, z) and the moments of the integrated intensity (Wj(z)} 
can be obtained from the normal characteristic function Cjy([3j,z) given in (14) in terms of the 
Laguerre polynomials for both single and compound (two-mode) fields (for details, see e.g. [|[ 
|l8| , [28| and references therein). Expressions for variances of quadrature components ((Apj(z)) 2 ) , 
((Aqj(z)) 2 ), for principal squeeze variance Xj(z), and for uncertainty product Uj(z) can also be 
found in [§ §§. 

We introduce only expressions for the variances ((AW) 2 } of the integrated intensity and for the 
principal squeeze variances A for single and compound modes, which will be used when analyzing 
the statistical properties on the basis of the short-length solution. 

We get the variance ((AW}) 2 ) of the integrated intensity in the form ([p|, p. 268) 

([AW,(z)] 2 > = B 2 (z) + |Q(z)| 2 + 2B j (z)Mz)\ 2 + {Cj{z)$\z) + c.c.) (17) 

for the jth single mode. The variance of the integrated intensity of the (j, fc)th compound mode 
equals 

<[AW, fe (z)] 2 ) = ([AW,(z)] 2 ) + ([AW k (z)] 2 ) + 2(AW ] (z)AW k (z)), (18) 
where the correlation of fluctuations of integrated intensities is given by 

(AW 3 (z)AW k (z)) = \D ]k (z)\ 2 + \D ]k (z)\ 2 

+ (D jk (z)C (z)C k (z) - D 3k (z)^(z)C k (z) + c.c.) . (19) 

The principal squeeze variances Xj(z) for the jth single mode and Xj k (z) for the (j, fc)th com- 
pound mode are expressed as J^] : 

\ 3 (z) = l + 2[B j (z)-\C j (z)\], 
\ ]k (z) = 2{l + B J (z) + B k (z)-(D jk (z) + c.c.) 

-\C 3 (z) + C k (z) + 2D 0k (z)\}. (20) 

We now can pay our attention to the statistical properties of the short-length solution. We 
firstly discuss the case with all modes being at the beginning of the interaction in coherent states 
(Brillouin processes). The coefficients of the normal characteristic function CV({/3j}, z) read up 
to z 2 : 

B Sl (z) = |<?sjV, 

B Vl (z) = \g Sl \ 2 z 2 , 

D Sl A 1 (z) = -g Al gs 1 z 2 /2, 
DstvAz) = igsi*, 

D Sl v 2 (z) = -gs 2 K* s z 2 /2. (21) 

The other coefficients (not obtained from the symmetry) are zero. We get from (17) and (21) for 
single mode variances of the integrated intensity 

«AW Sl (z)) 2 ) = 2|. 9Sl | 2 |£s 1 |V, 
((AW^(z)) 2 ) = 0, 

((AW^)) 2 ) = 2|. 9Sl | 2 |£v 1 |V. (22) 



7 



The symbols £ denote initial coherent amplitudes of the radiation and phonon modes. The variances 
of the integrated intensity for compound modes, which can be negative, read 

((AW Si aAz)) 2 ) = 2| ff5 J 2 |£ 5 jV-( 5Al5Sl ^ i G 1 * 2 + c.c.), 
((AVtWz)) 2 ) = mgs^kCv^ + c.c^ + lgs^il + ms^ + mvA 2 )^ 
+ (9s 1 9a 1 Ca 1 Cs 1 z 2 + cc.) + {g Sl ns€s 2 £v 1 z2 + c - c -)], 
((AW Sl v 2 (z)) 2 ) = 2| ff5l | 2 |C5 1 | 2 ^ + 2| ff 5 2 | 2 |^ 2 | 2 ^-(5s 2 4G i ey 2 ^ + cc.). 

(23) 

From the point of view of the short-length solution, negative variances of the integrated intensity 
reflecting nonclassical properties of light can be reached only in the compound modes (S\,Ai), 
(Si, Vi), and (Si, V2) when initial phases are suitably chosen. The linear coupling between Stokes 
modes (described by ks) leads to the negative intergrated intensity variance in the compound mode 
(Si, V2) and it also can support the negative integrated intensity variance of the mode (Si, Vi). 

The coefficients Cj (z) are zero not only in the short-length solution but also when we suppose 
initial non-squeezed light (Cj(0) = 0). This stemms from the momentum operator (1), which 
contains only first powers of operators. This excludes squeezing in single mode fields. Squeezing 
of vacuum fluctuations can be reached in the compound modes (Si, Ai), (Si, V\), and (Si, V2), the 
principal squeeze variances of which are given by, using (20) and (21): 

\siaM = 2[l + | ffSl |(| 3Sl |-k Ml |)z 2 ], 
\s lVl (z) = 2[\-2\g Sl \z +\g Sl \ 2 z 2 }, 

X Sl v 2 (z) = 2[1 + (\g Sl \ 2 + \gs 2 \ 2 )z 2 ~\gs 2 \\K S \z 2 }. (24) 

The linear Stokes coupling K5 thus can lead to generation of squeezed light in the mode (Si, V2). 

Now using the short-length solution we analyze statistical properties of modes for Raman 
processes, i.e. when phonons are initially chaotic and the other modes are initially coherent. 
Coefficients in the normal characteristic function Cj^({/3j}, z) have the form: 



B Sl (z) 


= l<7Si| 2 (^ +1)^ 2 , 


B Al (z) 




B Vl (z) 


= n Vl + |ffsi| 2 (nvi + l)^ 2 - \g A ^n v ^z 2 


D SiA 1 (z) 


= -9si9A 1 (nvi + f/2)^ 2 , 


DsivAz) 


= igsAnv! + l)^ , 


D S!V 2 (z) 


= -gs 2 K* s (n V2 + l)z 2 /2, 


Da 1 v 1 (z) 


= ig^nv^z , 


Da 1 v 2 (z) 


= 9*A 2 K An V2 z 2 /2. 



The other coefficients, which are not obtained from the symmetry, are zero. The mean numbers of 
initial chaotic phonons are denoted as n\r 1 and ny 2 for modes Vi and V2, respectively. 

Variances of the integrated intensity for both single and compound modes can be expressed 
using relations (17), (18) and (25): 

((A^ Sl (z)) 2 ) = 2\g Sl \ 2 (n Vl+ m Sl \ 2 z 2 , 

((AW Al (z)) 2 ) - 2\g Al \ 2 n Vl \UA 2 z 2 , 

((AW Vl (z)) 2 ) = n\ + 2\g Sl \ 2 n Vl (\£ Sl \ 2 + n Vl + l)z 2 

+ 2\g Al \ 2 nvMA 1 \ 2 -n Vl )z 2 + (2g* Ai g* Sl Zs 1 U 1 n Vl +c.c), 
((AW SlAl (z)) 2 ) = 2| 3Sl | 2 (n yi +l)|es 1 | 2 ^ 2 + 2|^ 1 | 2 ny 1 |CA 1 | 2 ^ 2 

- [9A 1 gs 1 (2n Vl + 1)£sAmZ 2 + c-c], 
((AW SlVl (z)) 2 ) = n 2 Vi +2\g Sl \ 2 (3\is 1 \ 2 + n Vl +l)(n Vl +l)z 2 
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+ [29A 1 9s 1 (n Vl + l^G^ + cc], 
((AW AlVl (z)) 2 ) = n 2 Vi +2\g Al \ 2 n Vl {nv 1 -\UA 2 )z 2 

- [2g Al g Sl n Vl C Sl C Al z 2 + c.c.]. (26) 

The above short-length expressions (26) indicate negative variances of the integrated intensity only 
in the compound mode (Si, Ai). The influence of linear Stokes (us) and anti-Stokes (k a ) couplings 
on the variances cannot be deduced from the expressions restricted up to z 2 . 

As has been discussed previously, squeezing of vacuum fluctuations cannot be obtained in 
single modes. However, the short-length solution provides interesting results for principal squeeze 
variances of the following compound modes, using (20) and (25): 

X SlAl (z) = 2[1 + \g Sl \ a (n Vl + l)z 2 + \g Al \ 2 n Vl z 2 - \g Al \\g Sl \(l + 2n Vl )z 2 }, 
*,SiVi(z) = 2[l + n Vl ~2\g Sl \(n Vl + l)z + 2\g Sl \ 2 (n Vl + l)z 2 - \g Al \ 2 n Vl z 2 }, 
X Al v 1 (z) = 2[l + n Vl - (ig* Al n Vl z + c.c.) + \g Sl \ 2 (n Vl +l)z 2 }, 
Xs 1 v 2 (z) = 2[1 + n V2 + \g Sl \ 2 (n Vl + l)z 2 + \gs 2 \ 2 (nv 2 + l)z 2 -\g A2 \ 2 n V2 z 2 

- \9s 2 \Ks\(nv 2 + l)z 2 ], 

X A t v 2 (z) = 2[1 + n V2 + \g s 2 \ 2 (n V2 + 1) z 2 + \g Al \ 2 n Vl z 2 -\g A2 \ 2 n V2 z 2 

- (g* A2 K A n V2 z 2 + c.c.)]. (27) 

Thus, squeezing can be reached in mode (Si,Ai) for small ny 1 - Especially, the above expresions 
show that Stokes linear coupling (ks) could generate squeezing in mode (Si, V2) and anti-Stokes 
linear coupling (k A ) could lead to squeezing in mode (Ai, V2). 

The advantage of the short-length solutions stemms from their simpler expressions which can 
be relatively easily analyzed. The relations above introduced enable us to reveal partly the role 
of Stokes and anti-Stokes coupling constants in generation of light with nonclassical properties. 
As it will be seen in the next section, the dynamics of the statistical properties is more complex 
in general than it may be seen in the short-length solution; but all tendencies are obvious from 
the short-length solution. The short-length solution is important especially for understanding the 
role of various initial phases, which can substantially change the statistical behaviour of nonlinear 
couplers. 



4 Statistical behaviour of nonlinear couplers 

This section contains a discussion of the statistical behaviour of the Raman and Brillouin couplers 
obtained on the basis of numerical resuls gained by applying the method described in section 2.1. 
In special cases, analytical solution from section 2.2 can be used to verify the correctness of the 
numerical approach. The discussion further generalizes results obtained in the previous section by 
means of the short-length solution. The couplers with Brillouin and Raman processes in operation 
will be discussed separately. 

Photon number distribution p(rij,z), reduced moments of the integrated intensity (W (z))/ 
(W(z)) k — 1, variances of quadrature components ((Apj(z)) 2 ), ((Aqj(z)) 2 ), principal squeeze vari- 
ance Xj(z), and the uncertainty product parameter Uj(z) (for definitions, see |gj) will be used 
when discussing the statistical properties of single and compound modes. Nonclassical behaviour 
of a field exhibits itself by sub-Poissonian photon number statistics, by negative reduced moments 
of the integrated intensity (they indicate antibunching of photons; a field with negative second 
reduced moment of the integrated intensity possesses a sub-Poissonian statistics) or by squeezing 
of quadrature variances (the variance is squeezed if its value is under 1 (2) for a single (compound) 
mode). 

Squeezing of vacuum fluctuations and sub-Poissonian statistics cannot be reached in single 
modes, if the fields are initially in coherent states or in coherent states with superimposed noise 
(this is connected with the fact that the coefficients Cj(z) in the normal characteristic function 
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Figure 2: Mean integrated intensity (W(z)) (solid curve without denotation) and reduced moments 
of the integrated intensity (W k (z)) /(W(z)) k - 1 for k = 2(*), k = 3(o), k = 4(A), and k = 5(o) 
for mode {S\,Ai)\ gg t — 1, gA x = 2, £si — 2i, £vi = 1, and the other parameters are zero. 



(14) are zero provided that they are initially zero). That is why we orientate our attention mainly 
to two-mode fields. 

We choose the condition \gA\ > \gs\ in cases discussed below, because such a condition is 
suitable for nonclassical light generation (19, 20 



4.1 Coupler with Brillouin processes 

We suppose initial coherent states in all six modes of the coupler. We divide the discussion to 
six sections; in each section a special configuration is discussed. These configurations change from 
simpler to more complex ones. 

Influence of k$ on stimulated Stokes process 

We suppose stimulated Stokes and spontaneous anti-Stokes processes in waveguide 1 (gs 1 7^ 0, 
(Mi 7^ 0, £,s! 7^ 0, £vi 7^ 0, £a x = 0). Reduced moments of the integrated intensity (in the following 
we simply speak of moments of intensity) for compound mode (Si , Vi ) can be negative if phases 
are suitably chosen (see the short-length solution (23)). Moments of intensity in modes (Si,Ai) 
and (Ai,Vx) are non- negative. Squeezing can be reached in modes (Si, A\) and (Si, Vi). There is 
a periodicity in the spatial development with the period 1 / \f\g~A 1 1 2 — \gsx I 2 (see [jl9| ) . Single and 
compound modes return periodically to coherent states (see Fig. 2 for mode (Si,Ai)). 

The introduction of the linear Stokes coupling (ks 7^ 0, £s 2 7^ 0) leads to the occurrence of 
negative moments of intensity also in modes (Si, Ai) (see Fig. 3a, compare it with Fig. 2), (S2, Vi), 
and (1S2, Ai). Fig. 3b illustrates a typical spatial development of the photon number distribution of 
compound modes with negative moments of intensity: the initially Poissonian statistics develope to 
regions exhibiting sub-Poissonian statistics in interchange with regions exhibiting super-Poissonian 
statistics. Squeezing of vacuum fluctuations can be reached also in modes (S2, Vi) and (S2,Ai) 
(see Fig. 3c for mode (S2, Ai)). The increase of ks decreases the period of the spatial dynamics 
(as indicated by the analytical solution in (11)) and it does not lead to an increase of noise for 
longer z, but to a tendency to conserve the initial statistics. 

The introduction of the Stokes process in waveguide 2 (gs 2 ^ 0, £,v 2 0) does not lead to 
nonclassical properties of mode (Vi, V2) at all. 

Influence of ka on stimulated anti-Stokes process 

We suppose stimulated anti-Stokes and spontaneous Stokes processes in waveguide 1 (gs 1 7^ 0, 
9 A x 7^ 0, £ytj 7^ 0, £vj 7^ 0, £s x = 0). Negative moments of intensity can be reached in mode 
(Si,Ai) if the initial phase of is suitably chosen (arg^J = ±n/2). Squeezing in modes 
(Si,Ai) and (Si,Vi) reflects nonclassical properties of fields in this case. Non-zero anti-Stokes 
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Figure 3: Reduced moments of the integrated intensity (W k (z)) / {W(z)) k — 1 for k = 2(*), k = 3(o), 
k = 4(A), and k = 5(o) for mode (Si, A\) (a), photon number distribution p(n, z) for mode (Si, Ai) 

(b) , and quadrature variances ((Ap(z)) 2 )(*), ((Aq(z)) 2 }(o), and principal squeeze variance X(z) (A) 

(c) for mode (S 2 ,Ai); g Sl = 1, fiUj = 2, k s = -10, £ Sl = 2, £ Vl = 1, £s 2 = 2 , and thc other 
parameters are zero. 
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Figure 4: The second reduced moment of the integrated intensity (W 2 (z)) / (W(z)) 2 — 1 (a) and the 
photon number distribution p[n, z) for mode {S2,A\) (b); gs 1 — 1, gA x = 2, ks = — 10, £si = 
^Ai = 2i, £vi = 1, £s 2 = 2, and the other parameters are zero. 



coupling (ka 7^ 0, £a 2 7^ 0) does not support nonclassical properties, but it can lead to transmission 
of light exhibiting negative moments of intensity from mode (Si, A\) to mode (Si, A2). Non-zero 
ka causes a fast increase of mean intensity and its moments. Faster oscillations in z occur in 
the spatial development of discussed quantities, but they are connected with a gradual increase of 
noise. There occur periods with a noise reduction. The coupling constant ka destroys gradually 
squeezing in modes (Si,Ai) and (Si, Vi). 

Influence of ks on stimulated Stokes and anti-Stokes processes 

Stimulated Stokes and anti-Stokes processes in waveguide 1 are assumed (gs t 7^ 0, gA t 7^ 0, 
^Si 7^ 0, 7^ 0, £vi 7^ 0). Negative moments of intensity and squeezing occur in both modes 
(Si, Ai) and (Si, Vi) if phases are suitably chosen (</>y +(f>s — 4>gs = — 7r /2, 4>v — 4>A + 4>gA — — 7r /2, 
see@). 

The interaction of mode S2 through Stokes coupling (ks 7^ 0, £s 2 7^ 0) leads to the occurrence 
of negative moments of the intensity and to squeezing also in modes (S2,A\) (see Figs. 4a and 
4b) and (S2,Vi). Fig. 4b demonstrates photon number distribution p(n,z) for mode (S2,Ai), 
which exhibits oscillations between sub-Poissonian and super-Poissonian statistics when z increases 
(compare Fig. 4b with Fig. 4a). 

If the phases are not suitably chosen, moments of mode (Si,j4i) are non-negative when the 
Stokes coupling is not included (ks = 0) (<fiy + (f>g — 4> gs = — 7r/2, (j>v — 4>a + <t>g A — tt/2, see Q). 
The inclusion of the Stokes coupling can stimulate negative moments of intensity in this mode. 

The addition of stimulated Stokes process in waveguide 2 (g$ 2 7^ 0, £s 2 7^ 0) does not lead to 
negative moments of intensity in mode (Vi, V2). 
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Influence of ka on stimulated Stokes and anti-Stokes processes 

Non-zero anti-Stokes coupling (ka 7^ 0, £a 2 7^ 0) of waveguide 2 suppresses nonclassical be- 
haviour of compound modes in waveguide 1 (gs 1 7^ 0, gA 1 7^ 0, £g t 0, £ai 0j ivi 7^ 0) for 
longer z. It introduces oscillations to the spatial development of quantities under investigation and 
it increases noise in moments of intensity and supports higher values of quadrature variances. How- 
ever, negative moments of intensity and squeezing for shorter z can be reached in mode (Si, A 2 ) 
as a consequence of the linear anti-Stokes coupling. 

The anti-Stokes process in waveguide 2 (gA 2 7^ 0, £y 2 7^ 0) does not create nonclassical effects 
in compound modes (Ai,A 2 ), (A 2 , Vi), (Ai,V 2 ), (Vi,V 2 ), and (Si,V 2 ). 

Influence of ks on stimulated Stokes and anti-Stokes processes in both waveguides 

Stimulated Stokes and anti-Stokes processes in both the waveguides are assumed now (g$ 1 =/= 0, 
9A l ^ 0, £ Sl ^ 0, ^ 0, Cvi ^ 0, g S2 ^ 0, g Aa ± 0, fa ? 0, U 2 ^ 0, £y 2 ? 0). Negative 
moments of intensity and squeezing occur in compound modes (Si,j4i), (Si,Vi), (S 2 ,A 2 ), and 
(S 2 ,V 2 ) when phases are suitably chosen. The linear Stokes coupling (k$ 7^ 0) conserves slightly 
negative moments of intensity and squeezing in these modes (see Fig. 5a for mode (Si,Ai)) and 
it can induce negative moments of intensity and squeezing also in compound modes composed of 
single modes in different waveguides, especially in modes (S 2 ,Ai), (S 2 ,Vi), (Si,A 2 ), and (Si, V 2 ) 
(see Fig. 5b for mode (S 2 , Vi); in mode (S%, A 2 ) substantial squeezing of vacuum fluctuations can be 
obtained) . The dynamics of the mean intensity and of moments in Fig. 5a shows a typical behaviour 
of compound modes; regions with slightly negative moments are followed by short regions of high 
increase of noise, when a decrease of the mean intensity changes to an increase. Complex values of 
Ks introduce asymmetry between the waveguides, which affects strongly the spatial development 
of the compound modes. 

Negative moments of intensity or squeezing are not generated in compound modes constituted 
by the same single modes in different waveguides, i.e. in modes (Si,S 2 ), (Ax,A 2 ), and (Vi,V 2 ) 
(see Fig. 5c for mode (Ai,A 2 )). However, these modes possess strong tendency to return to the 
coherent states. 

Influence of Stokes and anti-Stokes couplings on Brillouin processes in both waveguides 

Now we discuss a general configuration, i.e. we extend the configuration from the previous 
part by including anti-Stokes coupling (ka 7^ 0). Nonclassical effects in compound modes occur 
for the same modes, as discussed in the previous section, when values of the coupling constant ka 
were small. The increase of values of ka causes a successive increase of moments of intensity for 
both modes with negative moments (see Figs. 6a and 6b for mode (S 2 , Ai)) and modes with non- 
negative moments. This means that the occurrence of negative moments of intensity is restricted 
to shorter z for greater ka (compare Figs. 6a and 6b) or negative moments cannot occur at all. 

Greater values of ka mostly lead to a successive increase of values of variances and uncertainty 
product, but they can also serve to generate light, squeezing of which gradually develops with z 
(see Fig. 6c for mode (S\,A\)). 

If we compare the above results with the short-length ones, we can conclude that the general 
solution reveals richer possibilities for nonclassical light generation. The short-length solution 
together with the analytical one (l9) are extremely important for a suitable choice of phases, which 
are crucial for the generation of nonclassical light, as can be deduced from the short-length solution 
in the previous section. 

4.2 Coupler with Raman processes 

Initial coherent states in optical modes (Stokes, anti-Stokes) and initial chaotic states in phonon 
modes are assumed in the following discussion, which is again divided into the same six sections, as 
was done in the case of Brillouin processes. In general, the occurrence of regimes giving nonclassical 
effects in the statistical behaviour of light is similar as in the case of Brillouin processes. However, 
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Figure 5: Mean integrated intensity (W(z)) (solid curve without denotation) and reduced moments 
of the integrated intensity (W k {z)) / (W{z)) k - 1 for k = 2(*) and k = 3(o) for mode (Si, Ai) (a), 
reduced moments of the integrated intensity (W (z)} / (W {z)) k — 1 for k = 2(*) and k = 3(o) 
for mode (S2,Vi) (b), and mean integrated intensity (W(z)) (solid curve without denotation) 
and reduced moments of the integrated intensity (W k (z)} /(W(z)) k — 1 for k = 2(*), k = 3(o), 
k = 4(A), and k = 5(o) for mode (A 1 ,A 2 ) (c); g Sl = 1, gA x = 2, k s = 6i, gs 2 = 1, gA 2 = 2, 
^Si = ~ 2i, = 2i, £yi = 1, ^s 2 = — 2i, ^ 2 = 2i, £y 2 = 1, and the other parameters are zero. 



14 



(a) 



(b) 




(c) 



12 n 



10 



< 




Figure 6: Reduced moments of the integrated intensity (W k (z)) / (W(z)) k — 1 for k = 2(*) 
and k = 3(o) for mode (S2,A\) (a), the second reduced moment of the integrated intensity 
(W 2 (z)) / {W (z)) 2 — 1 for mode (S2, A\) (k a = i) (b), and quadrature variances ((Ap(z)) 2 )(*), 
((Aq(z)) 2 )(o), principal squeeze variance A(z)(A), and uncertainty product u(z)(o) for mode 
(Si, Ai) (c); ka — Qi and the other parameters are the same as in Fig. 5. 
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one can expect a diminished role of various phase relations for the occurrence of nonclassical 
light, owing to chaotic statistics of phonon modes, compared with the above discussed case. The 
short-length solution provides less information in this case. 

Influence of k s on stimulated Stokes process 

Stimulated Stokes process (gs x 7^ 0, qa x 7^ 0, £s x 7^ 0, n Vl 7^ 0) supports the squeezed light 
generation in compound modes (Si, ^i) and (Si, V\). Non-zero Stokes coupling (k s ^ 0, £s 2 0) 
supports nonclassical light generation. Especially, it induces antibunching in modes (Si,Ai), 
(Si,Vi), (S 2 ,Ai), and (S2,Vi) and also squeezing in modes (S2,^4i) and (S2,Vi). The increase 
of values of ny ± destroys antibunching in all modes. Antibunching in modes (S2, A\) and (S2, Vi) 
disappears for smaller ny 1 . 

Influence of ka on stimulated anti-Stokes process 

Squeezed light in compound modes (Si,Ai) and (Si, Vi) can be generated in stimulated anti- 
Stokes process (gs 1 7^ 0, gA x 7^ 0, £ s a 1 7^ 0, ny 1 7^ 0). If values of mean phonon number nv 1 are 
small, antibunching in mode (Si, A\) occurs. In this case, anti-Stokes coupling (ka 7^ 0, £a 2 7^ 0) 
leads to transmission of antibunching to mode (Si, ^2)- The anti-Stokes linear coupling supports 
also squeezing in mode (Si,A 2 ). However, antibunching and squeezing can occur only provided 
that n Vl is sufficiently small. 

Influence of ks on stimulated Stokes and anti-Stokes processes 

Negative moments of intensity in mode (Si, A\) and squeezing in modes (Si, A\) and (Si, Vi) 
characterize nonclassical properties of light for this case (gs 1 7^ 0, gA x 7^ 0, £si 7^ 0, 7^ 0, 
n v 1 7^ 0). Fig. 7a shows the photon number distribution p(n,z) for mode (Si,^4i), which de- 
veloped from Poissonian statistics through sub-Poissonian statistics to super-Poissonian statistics 
and this evolution repeats with increasing z. The corresponding moments of intensity, demon- 
strating photon antibunching at the corresponding z, are given in Fig. 7b. Squeezing of vacuum 
fluctuations in mode (Si, Vi) can be seen in Fig. 8. Increasing values of ny 1 smooth out negative 
moments of intensity. Thus, negative moments in mode (Si,Vi) cannot occur owing to chaotic 
phonon statistics (in the opposite to Brillouin scattering) and moments can evolve periodically and 
they can reach the initial values for longer z. But negative moments of intensity in mode (Si, Vi) 
can be induced by linear Stokes coupling (ks 7^ 0, £s 2 7^ 0). Non-zero Stokes coupling provides 
negative moments and squeezing in compound modes composed of single modes in different waveg- 
uides, especiallly in (S 2 ,Ai) and (S2,Vi). No antibunching and squeezing are possible in modes 
(Si, S 2 ) and (Vi, V 2 ) (gs 2 + 0, 9 a 2 + 0, n v . 2 ? 0). 

Influence of ka on stimulated Stokes and anti-Stokes processes 

Similarly as for Brillouin processes, non-zero anti-Stokes coupling (ka 7^ 0, £a 2 7^ 0) of waveg- 
uide 2 to waveguide 1 with active stimulated Stokes and anti-Stokes processes (gs 1 7^ 0, gA x 7^ 0, 
£si 7^ 0; 7^ 0, ny 1 7^ 0) results in gradual degradation of nonclassical properties of light, i.e. the 
corresponding values of moments of intensity as well as of the parameters characterizing squeez- 
ing successively increase. Non-zero anti-Stokes coupling also introduces oscillations to the spatial 
development of all the quantities under discussion, in agreement with the fact that the anti-Stokes 
interaction has tendency to conserve the initial photon statistics. 

Influence of ks on stimulated Stokes and anti-Stokes processes in both waveguides 

In this case, nonclassical effects occur in the same compound modes as for Brillouin processes 
discussed above, i.e. negative moments of intensity and squeezing can occur in compound modes 
(S 1 ,A 1 ), (Si,Vi), (S 2 ,A 2 ), (S 2 ,V 2 ), (S 2 ,Ai), (S 2 ,Vi), (S U A 2 ), and (S U V 2 ). 

Moments of intensity in mode (A\, V 2 ) cannot be negative, but values of moments for smaller 
z can be gradually reduced when z increases, as is shown in Fig. 9. 

Non-zero values of ks introduce oscillations to spatial development of discussed quantities. 
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Figure 7: Photon number distribution p(n, z) (a) and reduced moments of the integrated intensity 
(W k {z))/(W{z)) k - 1 for k = 2(*), k = 3(o), k = 4(A), and k = 5(o) (b) for mode 
= 1, = 2, = — 2i, = 2i, ny 1 = 0.1, and the other parameters are zero. 




Figure 8: Quadrature variances ((Ap(z)) 2 )(*), ((Aq(z)) 2 ) (o), and principal squeeze variance 
A(z)(A) for mode (Si, Vi); = 1 and the other parameters are the same as in Fig. 7. 
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Figure 9: Reduced moments of the integrated intensity (W k (z)) / (W (z)) k — 1 for k = 2(*), k = 3(o), 
fc = 4(A), and k = 5(o) for mode (A 1 ,V 2 ): g Sl = 1, fl 1 ^ = 2, k 5 = -6, gs 2 = 1, #a 2 = 2, £ Sl = -2i, 
= 2i, ny 1 — 0.1, £s 2 = 2, ny 2 = 0.1, and the other parameters are zero. 
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Figure 10: The second reduced moments (W 2 (z)) / (W (z)) 2 - 1 for modes (Si, V2) (*) and (£2, Vi) 
(°); 5Si = 1, 9a x = 2, k s = 6i, 3s 2 = !, 5a 2 = 2, = -2i, £ Al = 2i, n Vl = 0.1, £s 2 = -2i, 
= 2i, ny 2 =0.1, and the other parameters are zero. 

Complex values of k$ introduce asymmetry into spatial development of the statistical quantities, 
which affects mainly moments of intensity (see Fig. 10 for moments (Si,V^) and (S2,Vi)). The 
reduction of initial noise in compound modes can be obtained, including phonon modes V\ and V2. 

Influence of Stokes and anti-Stokes couplings on Raman processes in both waveguides 

The configuration from the previous section can be extended by including anti-Stokes coupling 
(ka 7^ 0). Negative moments of intensity and squeezing occur in the same compound modes as 
discussed in the previous section, provided that the values of the anti-Stokes coupling constant ha 
are small. Greater values of ka provide greater values of moments of intensity. Fig. 11 shows 
the spatial development of the second moment of intensity of mode (Si,A{), which is negative 
only for smaller values of z. Greater values of ha lead to a complete loss of negative values of 
moments for compound modes, in which they can occur for smaller values of ka- They also cause 
a successive increase of values of variances and uncertainty with increasing z in the most of modes, 
but squeezing for greater z can be also reached in some modes (e.g. (Si, A\), (S2, A2)). 

Greater values of ny 1 and ny 2 preserve negative moments of intensity and squeezing in com- 
pound modes excluding phonon modes ((S\,Ai), (Si,Ai), (S2,Ai), (S\,A 2 )). They destroy non- 
classical effects in compound modes involving phonon modes Vi and V 2 - 

Negative moments of intensity cannot be obtained in mode (Vi, V 2 ), but initial greater values 
of moments can be substantially reduced. 
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Figure 11: The second reduced moment of the integrated intensity (W 2 (z)) / (W (z)) 2 — 1 for mode 
(Si,Ai); ha — 6i and the other parameters are the same as in Fig. 10. 

If we compare the above results for couplers with Brillouin and Raman processes, we can con- 
clude that in both the cases sub-Poissonian photon statistics and squeezing of vacuum fluctuations 
can occur in both compound modes of one waveguide ((Sj, Aj), (Sj,Vj), j = 1,2) and compound 
modes consisting of single modes in different waveguides ((Si,A2), (Si, V2)). The Stokes linear 
coupling supports nonclassical properties of light, whereas the anti-Stokes linear coupling destroys 
them. Both Stokes and anti-Stokes linear coupling constants introduce oscillations to spatial de- 
velopment of the statistical quantities under discussion. 

5 Conclusions 

We have investigated nonlinear optical couplers composed of two waveguides operating by means 
of Brillouin or Raman processes, which mutually interact through linear Stokes and anti-Stokes 
couplings. We have solved the model, assuming classical strong laser pump modes, analytically 
for special cases and in short-length approximation as well. General solutions have been obtained 
numerically. We have investigated the statistical properties of single and compound modes, in 
particular we have obtained photon number distributions, reduced moments of the integrated 
intensity, quadrature variances, principal squeeze variance and uncertainty product. It has been 
shown that nonclassical light can be generated only in the compound modes. Light exhibiting 
sub-Poissonian photon statistics and squeezing of vacuum fluctuations can be obtained in both 
compound modes of one waveguide ((Sj,Aj), (Sj,Vj), j = 1,2) and compound modes composed 
of single modes in different waveguides ((Si, A2),(Si, V2)). The linear Stokes coupling supports 
nonclassical light generation, whereas the linear anti-Stokes coupling leads to its degradation. Both 
Stokes and anti-Stokes couplings introduce oscillations into spatial development of the statistical 
quantities. The above conclusions are valid for both Brillouin and Raman processes. However, 
in Raman processes, increasing values of the mean phonon numbers degrade nonclassical light, 
especially in compound modes involving phonon modes. 
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